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ABSTRACT 

X> ■ 
<D ■ 

pL| ' Previously, we showed that, owing to effects arising from quantum electrodynamics (QED), 

£ — ■ magnetohydrodynamic fast modes of sufficient strength will break down to form electron-positron 

i— I \ pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of 

the energy of the fast mode fuels the development of an electron-positron fireball. However, 
j> • a small, but potentially observable, fraction of the energy (~ 10 33 ergs) can generate a non- 

\ thermal distribution of electrons and positrons far from the star. In this paper, we examine the 

cooling and radiative output of these particles. We also investigate the properties of non-thermal 

■ emission in the absence of a fireball to understand the breakdown of fast modes that do not yield 
an optically thick pair plasma. This quiescent, non-thermal radiation associated with fast mode 
breakdown may account for the recently observed non-thermal emission from several anomalous 
X-ray pulsars and soft-gamma repeaters. 

1. Introduction 

£3 ■ 

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) comprise a subclass of neutron stars 
with unusual properties compared to typical radio pulsars. The dozen or so known SGRs and AXPs emit 
pulsed X-rays, with steadily increasing periods, as well as bursts of hard X-rays and soft 7-rays at irregular 
intervals, spin slowly, and appear to be radio quiet. Compared with high-mass X-ray binaries, SGRs and 
AXPs have relatively low X-ray luminosities, soft spectra, and no detectable companions (for reviews, see, 
e.g. Hurley 2000; Mereghetti et al. 2002). 

While it seems clear that SGRs and AXPs cannot be rotation-powered, their ultimate energy source 
remains unclear. In the "magnetar" model (Duncan & Thompson 1992; Thompson & Duncan 1995), these 
objects are ultramagnetised neutron stars with fields £?* ~ 10 14 — 10 15 G, fueled by either magnetic field 
decay (Thompson & Duncan 1996; Heyl & Kulkarni 1998) or residual thermal energy (Heyl & Hcrnquist 
1997c, d). Alternatively, these stars may have "ordinary" field strengths £?* ~ 10 11 — 10 13 G and be powered 
by accretion from e.g. a fossil disk (Corbet et al. 1995; Chatterjee et al. 2000; Alpar 2001; Menou et al. 
2001), but recent analyses indicate that disks of the required mass may not survive sufficiently long around 
these stars (e.g. Ek§i et al. 2005). 

According to the magnetar model, bursts occur when fractures in the crust of the star send Alfven 
waves into the magnetosphere (Thompson & Duncan 1995). In earlier work, we analysed wave propagation 
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through fields exceeding the quantum critical value -Bqed = m 2 c 3 /eh 4.4 x 10 13 G, and demonstrated 
circumstances under which electromagnetic (Heyl & Hernquist 1998) and some MHD waves, particularly fast 
modes (Heyl & Hernquist 1999) evolve in a non-linear manner and eventually exhibit discontinuities similar 
to hydrodynamic shocks, owing to vacuum polarisation from quantum electrodynamics (QED). In Heyl & 
Hernquist (2005, hereafter Paper I), we developed a theory to account for bursts from SGRs and AXPs based 
on "fast-mode breakdown," in which wave energy is dissipated into electron-positron pairs when the scale 
of these discontinuities becomes comparable to an electron Compton wavelength. We showed that, under 
appropriate conditions, an extended, optically thick pair-plasma fireball would result, radiating primarily in 
hard X-rays and soft 7-rays. 

Our mechanism provides a natural and efficient outlet for ultramagnetised neutron stars to form bursts 
of high energy radiation without significant emission at shorter wavelengths. However, as we noted in Paper 
I, non-thermal electrons and positrons would be produced far from the star, beyond the optically thick 
fireball, radiating over a wide range of frequencies. In addition, MHD fast modes of insufficient amplitude to 
generate an optically thick fireball will still dissipate through pair-production, seeding non-thermal emission. 

In what follows, we extend our previous investigation of fast-mode breakdown to estimate the spectrum 
of non-thermal emission expected outside the region containing an optically thick fireball. We also consider 
the fate of fast modes which dissipate their energy through pair production but at a rate insufficient to 
yield a fireball. We show that when a fireball is generated with sufficient energy to account for the thermal 
radiation associated with bursts from SGRs or AXPs, small, but detectable amounts of non-thermal emission 
are predicted which dominate the spectrum both well below and well above X-ray energies. Our model makes 
a specific prediction for the spectral energy distribution of the non-thermal component that can be used to 
test our theory for SGR and AXP bursts. 

For weaker fast modes that do not yield an optically thick fireball, we show that the non-thermal 
radiation may be sufficient to explain the quiescent, non-thermal emission observed recently from several 
SGRs and AXPs. Therefore, our model naturally provides a unified explanation for several of the unique 
characteristics of SGRs and AXPs in the context of the magnetar scenario. 

2. Calculations 

In Paper I, we assumed that the energy dissipated from the fast mode fuels a pair fireball. This 
assumption is reasonable as long as the region is opaque to X-rays. Any electrons and positrons produced 
will generate synchrotron radiation which will be scattered and form additional pairs in the strong magnetic 
field. Beyond about one hundred stellar radii, the synchrotron radiation can escape without becoming 
thcrmaliscd. Although it may suffer synchrotron self-absorption, a power-law distribution of electrons and 
positrons can radiate outside the fireball, and we would expect a synchrotron radiation spectrum. 

2.1. Primary Pairs and Photons 

Below, we examine our model from Paper I in the limit when the pair plasma does not thermalise. 
As fast modes produced near the stellar surface propagate outwards, they develop strong discontinuities 
(shocks) when the magnetic field is sufficiently large and the wave energy will eventually be dissipated 
through pair production. If we consider the situation when the electron-positron distribution is no longer 
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thermal, we must adopt some model for the energies of these particles when they are created. A simple 
and reasonable assumption is that they are nearly at rest in the frame of the shock as they are produced. 
In the ultrarelativistic limit, this is required in order to conserve both momentum and energy as the wave 
dissipates. 

The shock is ultrarelativistic and travels at the speed of light in the strong magnetic field. The non- 
thermal pairs are produced far from the star where it is appropriate to use the weak-field limit of the index 
of refraction, 

where a/ = e 2 /ch is the fine-structure constant, and so 

- 1 _ 1 _ 3 / 57r -Bqed 

where Sqed = m 2 c 3 /(eh) w 4.4 x 10 13 G. (We note that the derivation of the effective Lagrangian of QED 
presented by Heyl & Hernquist (1997a), valid for any magnetic field strength, correctly reduces to the weak 
field limit of e.g. Heisenberg & Euler (1936), yielding the above expression for the index of refraction (Heyl 
& Hernquist 1997b).) 

Paper I describes in detail how the energy is dissipated as the fast mode travels away from the stellar 
surface. Here, we are interested in the evolution of the wave in the weak-field limit, far from the star, so 
we can obtain simple relationships for the amount of energy released as a function of the distance from the 
surface. The emission from outside a particular radius is proportional to (kb)~ 2 r~ 3 where k is the initial 
wavenumber of a fast-mode wave, b is its initial amplitude compared to the magnetic field of the star, and 
r is the distance from the star, so we have 

dE = Ar~ A dr , (3) 

where the constant A must be determined from modeling the evolution of the wave near the neutron star, 
as described in Paper I. 

At each radius from the star, we assume that the pair spectrum is monoenergetic, so the total emission is 
the sum of contributions from different electron-positron energies traveling in different strengths of magnetic 
field. The particles more distant from the star have higher typically energies and also travel through weaker 
fields. 

We explicitly add the radial dependence to the Lorentz factor of the particles and calculate the Larmor 
radius of the particles as a function of distance from the star. 

where R is the radius of the star, and B* is the pole field strength at the stellar surface, so 
-fmc 2 3 /57T ( B \~ 2 h oe m-15 v>-t( r \ h ( ^ 



farmer « = ~4 - - — = 2.6 X lO" 10 ^, 
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where i? 6 is the radius of the star in units of 10 6 cm, and we see that particles produced closer than a 
thousand stellar radii are essentially trapped. From Paper I, we find that the fireball extends to about one 
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hundred stellar radii, so our region of interest spans an order of magnitude in distance but three orders of 
magnitude in magnetic field strength and particle energy. 

If classical synchrotron radiation is important, the emission at each radius will extend approximately 
up to photon energies of 

fc 2 fc 2 eB 45 7T i?QED 2 fa\ 

hw = "ffrw B =T — ~ -; tt - mc , (6) 

mc 4 a/ B 

where Hujb is the non-relativistic cyclotron energy. 

We can estimate the recoil on the electron owing to the emission of individual photons 
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and, so, we must treat the synchrotron emission quantum mechanically. 

According to the results of Bcrcstctskii ct al. (1982), the emission peaks at a energy 

77 1 fe 1 2 X 1 2 3 5w B QED 2 

2 2 |+ x 2 4\j a f B 

The power radiated by the electron is given by 

32r(|) / x2/3 (mc 2 ) 2 (mc 2 ) 2 

so the electrons and positrons lose their initial energy over ~ 2O7 Compton wavelengths and the primary 
radiation is strongly beamed parallel to the direction of the fast wave. 



2.2. Formation of Secondary Particles 

The high-energy synchrotron photons passing through the strong magnetic field may produce pairs 
through one-photon pair-production. It is reasonable to assume that Hco ;§> mc 2 (Eq. 8), so we will use the 
results of Berestetskii ct al. (1982). The key parameter is essentially 

H 2 eBw B huj x 3 /5~t 

k= — 5-r- = — 5- w - = -\ — w 35 > 1 , (10) 

m 6 & Bqed mc^ 2 4y a/ 

where we have used the value of ui c from Eq. (6). Therefore, we can use the large k limit of the pair-production 

_ 5r 2 (|) 3 Va g mc 2 3 10 W^ (n) 

These secondary pairs will typically have energies one-quarter that of the first generation of particles. The 
fourth generation of particles will have energies typically 4 -3 or 64 times smaller than the first generation. 
This fourth generation will have \ ~ 1 so the typical energy of the resulting synchrotron photons will be less 
than one- half the energy of the pairs so k <C 1, quenching the pair production rate, 

3 3 / 2 B mc 2 _ S/3K 

w = ¥T^B^^ e 7 (12) 
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We assume that that c/w <r for pair production to be effective. We find that the pair-production quenches 
for 

r B 



K ~ O.f 



1 + 0.04 In 



(13) 



^lO 8 cm 10- 5 B Q e D/ 

To calculate the final spectrum of radiation, we will make the ansatz that only the final generation of pairs 
contributes to the observed emission and that their energies are about half the energy of the last generation 
of pair-producing photons, so they have \ w n/2 ~ 0.04. Using the definition of \ we nn( l 



7 = X - = 0-05- v 



B 



B 



(14) 



If we compare this result with Eq. (2), we find that the energy of the pairs has been degraded by a factor of 
about 1700. 



Because \ ^ 1 we can assume that the emission is classical. The emission will have typical energies of 

hw t = ^hw B «2.5 x 10- 3 ^^mc 2 . (15) 



We can combine Eq. (3) with Eq. (15) to obtain an estimate of the spectrum of resulting photons, 

dE oc r- 6 dE 7 oc E~ 2 dE 1 for £ 7 > £ brC ak . (16) 
where -Ebroak is the value of Eq. (15) at the edge of the fireball. 

A final important element is the cooling time of the final generation of pairs: 



jmcr 
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where r = e 2 /mc 2 is the classical electron radius, and so the final generation of pairs cools nearly immedi- 
ately The pairs closest to the surface of the star have the lowest energy and cool the most quickly, so for 
lower energies we have 

dE oc E- 1 ' 2 dE 1 for E < E 1 < £ broak , (18) 

where Eq = kujB,FB is the non-relativistic cyclotron energy at the outer edge of the fireball. Below this 
energy we observe radiation from the pairs formed further away from the surface of the star. There are fewer 
of these pairs, but because they form in a weaker magnetic field their spectra will extend to lower energies. 

The spectra of these electrons extends as E~ x l 2 from i?break Bfb/B down to EqB/Bfb but the normal- 
isation is smaller by (B / B^b) 2 so we have 



dE oc E 1 dE^ 1 for E n 



<E 7 <E , (19) 
where i? m in is the cyclotron energy at the radius where the pairs are no longer trapped from Eq. 5. 



Putting together Eq. (16), (18) and (19) and normalising by the total energy in non-thermal radiation 
yields the complete spectrum 
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(20) 
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where E m i n and E maK are determined by the non-relativistic cyclotron frequency and the typical emission 
energy, given in Eq. 16 at the outer edge of the region where the primary pairs are initially trapped by the 
field from Eq. 5. Typically, E min ~ 10~ 3 E and -E max ~ 10 3 i?broak- To find the normalisation in Eq. (20), 
we have assumed that E max 3> -EWcak and E -C inbreak 

The bulk of the energy emerges near -Etreak, and the bulk of the photons have energies near Eq, and 
the mean energy of a photon is 

V^brcak^o = xmc 2 = 0.05mc 2 , (21) 

where \ is given in Eq. (14). 

We now estimate the values of Eq and -EWeak- From the models presented in Paper I, we have rpe « 
100R, so 

B B 
E « lO" 6 — mc 2 w 5 eV and £ bre ak ~ 1600— ^-mc 2 w 80 MeV , (22) 

-DQED B* 

where the final results assume that £?* = 10-Bqed- The mean photon energy is around 25 keV. In Paper I, 
we argued that about 10 ergs per burst may go into non-thermal emission. In this case, we would have 
approximately 2.5 x 10 41 photons or 2 x 10~ 3 photons per square centimeter at a distance of 10 kpc. The 
bulk of the photons will have energies around E . 

The total energy in the non-thermal emission is proportional to the value of the magnetic field at the 
edge of the fireball (from Eq. 3), so as one can see from Fig. 1, the high-energy spectrum simply extends 
to a lower value of -Ebreak and the low-energy spectrum extends to a higher value of Eq as the radius of the 
fireball decreases. 



2.3. Time Dependence of Spectral Emission 

Fig. 1 gives the total spectrum of radiation, both thermal and non-thermal from the burst, integrated 
over time. While the pairs are still relativistic we find that the cooling time is 

3 / E B \ J% 

Tcool = ~ o ~d 2 ( 23 ) 

2 \mc 2 Bq EB J mc 2 

as a function of the energy of the typical synchrotron photon E t . If we compare this result with Eq. (15), 
we find that initially the typical photon energy is also inversely proportional to the field strength, so as time 
passes we have a spectrum Fe oc E~ 2 above a frequency E^ ca ^(t) oc t~ 2 . The value of Fe at -Ebreak(i) is 
independent of time. 

The pairs produced nearest to the star become non-relativistic first so as inbreak approaches Eq the 
spectrum steepens and at late times approaches Fe oc _E 7 , but now the normalisation decreases rapidly with 
time. To appreciate the rate of cooling we can substitute a value for Eq into the cooling time Eq. (23) to get 

-co ol = !^A (24) 
2 B mc 1 

so the pairs cool very rapidly to non-relativistic energies over 10 5-8 Compton times or 0.1 to 10 femto- 
seconds. Observationally, one expects to see the non-thermal spectrum as depicted in Fig. 1 immediately 
before the fireball radiation (~ 10~ 3 s). The non-thermal emission dominates the fireball emission both at 
low energies and for E 1 ^> kTpB 
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5 
log 10 E [keV] 

Fig. 1. — Spectrum from an SGR or AXP fireball (smooth solid curve) compared with the non-thermal 
emission. Tfb — 10 keV and the non-thermal emission comprises 10 -5 of the fireball energy. The value 
of the magnetic field at the edge of the fireball (100 stellar radii) is taken to be 3 x 10~ 5 £?qed, yielding 
Eq = 15 cV and -Ebrcak = 40 MeV (broken solid line). If the fireball ends at thirty stellar radii, Eq and 
-Ebreak increase and decrease respectively by a factor of 25000 and the total non-thermal energy is 3 x 10~ 4 
of the fireball emission (dashed line). As we argue in § 2.7, values of inbreak less than 1 MeV do not make 
sense physically, so the fast- mode cascade as we have described is limited to fields less than « 1CT 3 -Bqedj 
well within the weak- field limit. 
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2.4. Inverse-Compton versus Synchrotron Power 

We have assumed implicitly that synchrotron emission dominates the energy loss of the high energy 
pairs produced in the cascade. If we assume that the fireball emits 10 40 erg/s we can estimate the photon 
energy density: 

photon - -^r c = 2 x 10 12 L 40 r 8 - 2 erg cm" 3 (25) 



and the magnetic energy density 



so 



U B = ^ = Ax 10 16 B 2 15 r 8 - 6 erg cm" 3 (26) 

= t^ton = g x 1Q -6 i4or 4 B 2 5 (2?) 
-TSych L'B 

Inverse-Compton emission may be important for r 3> 10 9 cm. Coincidcntally, this is the same distance from 
the star where the Larmor radius of the initial pairs becomes greater than the distance from the star (Eq. 5), 
where synchrotron emission is not very effective. Regardless, the total energy deposited in this outermost 
region is 10 -3 lower than the synchrotron-dominated region (Eq. 3). 

In principle, the larger bursts such as the March 7 and August 28 events have sufficient luminosity to 
make inverse Compton emission important. However, the cooling time for the pairs (Eq. 17) is so short that 
the pairs will lose most of their initial energy before the fireball radiation reaches the non-thermal emitting 
region. The cooling time for pairs emitting at lower energies is sufficiently long that inverse- Compton emission 
might be important for these electrons in the brightest SGR bursts. 



2.5. Two-photon Pair Production 

To determine the efficacy of two-photon pair production we must first decide whether the highest energy 
photons produced in the cascade (Eq. 8) lie above the threshold to interact with the thermal photons from 
the neutron star or other ambient photons. The photons produce by the fireball simply do not catch up with 
the photons produced in the cascade. 

For the thermal photons, we must have 

E CM = y^pcakffthcrmai (1 - cos 9) « 0.133 ^g™ ' MeV > 1 MeV (28) 

so for 

r > 60-^- ( £th ; mal ) R (29) 
-Bqed V 1 keV / 

two-photon pair production may operate. However, the cascade typically dumps most of its energy around 
r = lOOi? and £?* > 10-Bqed so two-photon pair production against the thermal photons from the neutron 
star lies below threshold. Furthermore, even if the thermal photons were above threshold, the mean-free 
path is greater than 10 3 R and much much greater than the mean-free path for one-photon pair production. 

We can repeat the calculation for ambient photons from starlight 

ECU = ViWSambicnt = 4 x I0' 3 ^ ^) ^ MeV (30) 
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so ^ 

r>4 o(-^) p ambicnt N l R. (31) 
\BqedJ \ 1 cV J 

The ambient starlight typically does have sufficient energy to pair-produce against the first generation of 
photons produced in the cascade but the mean-free path is from kiloparsecs to megaparsecs much greater 
than the scale of the cascade region. 



2.6. The Fast-Mode Cascade Contrasted with the Pulsar Cascade 

The main difference between the fast-mode cascade and the pair-cascade invoked to explain pulsar 
emission is geometry. The pulsar cascade is instigated by charged particles that are restricted to travel along 
magnetic field lines - any perpendicular momentum is quickly reduced through synchrotron emission. The 
photons that are produced also travel nearly parallel to the field lines, so photons must travel some fraction 
of the curvature length of the magnetic field before their momentum perpendicular to the magnetic field is 
sufficient to produce a pair. Consequently, the pairs are produced nearly on threshold. 

The fast-mode is not restricted to travel along the field lines. Here, we have explicitly assumed that 
it travels perpendicular to the field lines. Regardless of the precise angle, the pairs produced initially will 
have much of their momentum perpendicular to the field resulting in many photons traveling perpendicular 
to the field that can in principle produce pairs in situ. What quenches the process is the reduction of the 
cross-section as the energy of the subsequent photon generations decreases. 



2.7. The Fast-Mode Cascade without SGR Bursts 

In Paper I, we argued that a sufficiently widespread production of fast modes would result in the creation 
of a fireball because the vicinity of the star would be opaque to X-rays simply owing to the Thomson cross 
section of the pairs. The optical depth between the surface of the star and infinity is simply proportional to 
the fraction of the star's surface area generating magnetohydrodynamic waves. The SGR bursts typically 
involve a large surface area ~ R 2 emitting fast-modes with a fractional amplitude of several percent within 
several light-crossing times of the inner magnetosphere. 

These large-scale events are bound to be rare. On the other hand, smaller scale stresses in the crust of 
the star may be relieved more regularly. This would excite only a small fraction of the neutron star surface, 
but in principle the amplitude of the waves might be sufficient to cause fast-mode breakdown. In this case, 
all of the radiation produced would follow the non-thermal distribution outlined in this paper. 

What would be the typical flux of these small-scale fast modes? Thompson & Duncan (1996) and 
Heyl & Kulkarni (1998) have argued that the quiescent emission of SGR and AXP neutron stars may be 
powered by the decay of the magnetic field. The quiescent thermal emission may only be a small fraction 
of the total energy released by the decay of the magnetic field. Recent observations of SGRs and AXPs 
indicate that the thermal radiation may indeed be just the tip of the iceberg (Molkov et al. 2004; Mereghctti 
ct al. 2004a; Kuiper et al. 2004); therefore, we will be liberal in our assumption of the relative amplitude of 
the thermal and non-thermal radiation and allow the pair cascade to operate beyond a certain radius from 
the star or equivalently below a certain magnetic field strength (-B max ) • This model has two parameters, 
-B max determines the position of the two breaks in the spectrum and the total normalisation. The process 
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of fast-mode breakdown predicts a particular relationship between the location of the two breaks in the 
spectrum (22) and particular slopes (20) between and beyond the breaks. 

Figure 2 shows the observed broad-band spectrum of several AXPs and SGRs. To compare the spectra 
we have placed the more distant objects (AXP IE 1841-045 and SGR 1806-20) whose hard X-ray emission 
was discovered with INTEGRAL at the distance of 4U 0142+61 whose optical emission (Hulleman et al. 
2000) is very likely to be nonthermal (Ozel 2004). 

Because the location of the two breaks in the spectrum both depend on the strength of the magnetic 
field at the inner edge of the breakdown region, we find that the presence of extensive nonthermal optical 
emission indicates that the non-thermal hard X-ray emission should peak at about 30 MeV, a factor of two 
hundred beyond the observed spectrum. The best limits in this energy range are provided by EGRET. 

We come to this conclusion under the assumption that SGR 1806-20 and AXP IE 1841-045 have a similar 
optical excess to 4U 0141+61. A more conservative assumption would be that the hard X-ray emission does 
not extend far beyond the observations from INTEGRAL with spectral breaks at about 1 MeV and 650 eV. 
This situation is somewhat natural. The fast-mode cascade is limited to pairs with sufficient energy to 
produce photons with E > 1 MeV that can subsequently pair produce. Lower energy electrons simply cool, 
giving the observed cooling spectrum in the hard X-rays. The total energy in the non-thermal emission is 
also reduced by a factor of a few. In the context of the fast-mode cascade it is difficult to have inbreak < 2mc 2 . 
This model is denoted as the "Minimal Model" in Fig. 2 because -Ebreak takes on the minimal value that 
makes sense physically; i.e w 1 MeV. 

We can also obtain a conservative picture for AXP 4U 0142+61 by accounting for only the optical excess 
and getting a prediction for the hard X-ray emission. We assume that E is a few electron volts giving inbreak 
on the order of 200 MeV. The minimal hard X-ray emission required to account for the optical emission is 
depicted in Figure 2 as the "Optical Model" . This optical model agrees well with the results reported by den 
Hartog et al. (2004). On the other hand, if AXP 4U 0142+61 has inherent non-thermal emission as strong 
as either AXP IE 1841-045 or SGR 1806-20, the EGRET flux would be larger as given under the "Unified 
Model" in Table 1. 

Table 1 gives the predicted EGRET fluxes above 100 MeV for the various objects and models. The 
optical model and the unified model predict approximately similar EGRET fluxes for 4U 0142+61. The 
minimal model based solely on the INTEGRAL data exhibits a flux above 100 MeV about two hundred 
times smaller than the unified model. Because the optical model cannot explain the observed INTEGRAL 
data for IE 1841-045 and SGR 1806-20, it is omitted. Similarly, the minimal model cannot explain the 
optical data for 4U 0142+61. We sec that for IE 1841-045 and SGR 1806-20 the predictions for the minimal 

Table 1: Predicted flux above 100 MeV and observed EGRET upper limits (Grenier & Perrot 2005) in units 
of 10~ 8 photons s _1 cm -2 . The GLAST upper limits are nominally (0.2 - 0.4) x 10~ 8 photons s _1 cm -2 
(GLAST Facilities Science Team 1999). 
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Fig. 2. — The spectrum produced by fast-mode breakdown is superimposed over the observed thermal and 
non-thermal emission from several AXPs and SGRs for models that fit either the optical or INTEGRAL 
data solely and one that fits both sets of data. The unabsorbed optical data are from Hulleman et al. (2000) 
via Ozel (2004) for AXP 4U 0142+61. The uppermost black symbols are the hard X-ray band are from 
Molkov et al. (2004) for SGR 1806-20. Mereghetti et al. (2004a) obtained similar results for the SGR. The 
middle sets of points in the hard X-ray data (blue is total flux and red is pulsed flux) are from Kuiper et al. 
(2004) for AXP IE 1841-045. The green squares plot the INTEGRAL data reported by den Hartog et al. 
(2004) for AXP 4U 0142+61. We normalised the den Hartog et al. (2004) results using the observations of 
the Crab by Jung (1989). We scaled the emission from the three sources by assuming that they all lie at 
the distance of AXP 4U 0142+61. We used 3 kpc for AXP 4U 0142+61 (Hulleman et al. 2000), 7.5 kpc for 
AXP IE 1841-045 (Sanbonmatsu & Helfand 1992) and 15 kpc for SGR 1806-20 (Molkov et al. 2004). 
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model lie comfortably below the EGRET upper limits. In the context of the fast-mode breakdown model, 
this means that the optical emission for IE 1841-045 and SGR 1806-20 is inherently weaker than from 
4U 0142+61. 

On the other hand, 4U 0142+61 is difficult to explain in the context of cither model because of its 
large optical flux. Perhaps 4U 0142+61 was more active during the epoch of the optical observations than 
during the EGRET observations. Some AXPs exhibit variable X-ray emission such as AX J1845-0258 and 
IE 1048.1-5937 (Vasisht et al. 2000; Mereghetti et al. 2004b) so this conclusion might be natural. 

3. Discussion 

We have presented a model for the non-thermal emission from AXPs and SGRs. It is a natural extension 
of the model for bursts that we presented in Paper I. When a dislocation of the surface of the neutron star 
is sufficiently large, the resulting fast modes will produce sufficient pairs to make the inner magnetosphcrc 
of the neutron star opaque to x-rays, generating a fireball. Some small but observable fraction of the 
energy initially in the fast modes is dissipated outside the opaque region yielding a characteristic fast-mode 
breakdown spectrum. A second possibility is that the crust of the neutron star is constantly shiftly over 
small scales, generating fast modes whose breakdown is insufficient to produce a fireball. In this case we 
would associate the non-thermal radiation with the quiescent thermal radiation from the surface of the star. 

The exponential quenching of one-photon pair production and the dipole geometry of the neutron star 
field at large distances constrains the spectrum produced by fast-mode breakdown to have only two free 
parameters. The cyclotron energy (E ) at the inner edge of the breakdown region (or alternatively £"break) 
and the normalisation of the spectrum. The model does not yield any freedom in setting the spectral slope. 
Below E , EFe is proportional to E. Between E and .Ebreak it is proportional to E 1 / 2 and above -Ebroak 
EF E oc E' 1 . 

Because of the paucity of free parameters, our model has great predictive power, if one can ascribe the 
observed flux from an AXP or SGR to fast-mode breakdown. In the case of 4U 0141+61 whose optical and 
near infrared flux has been determined (Hullcman et al. 2000) , the hard X-ray flux from fast-mode breakdown 
should be at least as large as found by den Hartog et al. (2004). Additionally, we have several predictions 
for the gamma-ray flux from AXPs and SGRs. Regardless of which model one assumes, the emission should 
extend beyond observed INTEGRAL data without a break below 1 MeV. Further INTEGRAL observations 
of these objects may verify this claim. If the object has a significant optical excess, such as 4U 0141+61, 
we predict that EFe should continue to rise to 10 — 200 MeV, well into the realm of GLAST. Even the 
least favourable "Minimal" model for AXP IE 1841-045 and SGR 1806-20 yields a marginal detection with 
GLAST. 

Observations of the AXPs and SGRs continue to surprise, as do theoretical investigations of ultramag- 
netised neutron stars. In this paper and Paper I we have presented a unified model for the thermal burst 
emission and non-thermal emission from ultramagnetised neutron stars. The model has few underlying as- 
sumptions: magnetars produce fast modes sufficient to power the non-thermal emission and, more rarely, 
the bursts, the magnetic field far from from the star is approximately dipolar and quantum electrodynamics 
can account for the dynamics of pairs and photons in strong magnetic fields. The model for the non-thermal 
emission has two free parameters, a normalisation and break frequency. Further observations can easily 
verify or falsify this model and potentially provide direct evidence for the ultramagnetised neutron stars that 
power AXPs and SGRs and the macroscopic manifestations of QED processes that account for their unique 
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attributes. 

J.S.H. would like to thank Alice Harding for help with the EGRET upper limits. J.S.H. is supported 
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